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We study the relation of confinement and chiral symmetry breaking in gauge theories with non- 
trivial center, such as SU(iV) gauge theories. To this end, we deform these gauge theories by 
introducing an additional control parameter into the theory and by varying the representation of 
the quark fields. We then consider a large-d(R) expansion of the effective action, where d(R) denotes 
the dimension of the representation R of the quark fields. We show how our large-d(R) expansion can 
be extended in a systematic fashion and discuss the effects of l/d(R)-corrections on the dynamics 
close to the finite-temperature phase boundary. Our analysis of the fixed-point structure of the 
theory suggests that the order, in which the chiral and the deconfinement phase transition occur, is 
dictated by the representation of the quark fields and by the underlying gauge group. In particular, 
we find that the phase diagram in the plane spanned by the temperature and our additional control 
parameter exhibits an intriguing phase structure for quarks in the fundamental representation. For 
SU(iV) gauge theories with adjoint quarks, on the other hand, the structure of this phase diagram 
appears to be less rich, at least in leading order in the I/d(R)-expansion. 



I. INTRODUCTION 

For quantum chromodynamics (QCD) with iV = 3 col- 
ors and two (light) quark flavors it has been found in vari- 
ous studies that the phase transitions associated with chi- 
ral symmetry restoration and deconfinement lie remark- 
ably close to each other, see e. g. Refs. [HUB], From 
a phenomenological point of view, this observation has 
important consequences for our understanding of the dy- 
namics in heavy-ion collisions as well as of the generation 
of hadron masses in the early universe. In fact, a com- 
prehensive picture of the dynamics close to the finite- 
temperature phase boundary of QCD is required for a 
reliable description of data from heavy-ion collision ex- 
periments [T7] , 

While the underlying mechanisms associated with the 
confinement of quarks are not yet fully understood, we 
have a profound understanding of chiral symmetry break- 
ing in gauge theories. In fact, it is already known from 
Nambu-Jona-Lasinio (NJL) models [HI Qjj] that chiral 
symmetry breaking is indicated by strong quark self- 
interactions. In particular, the four-quark interactions 
play a prominent role since they are directly related 
to the chiral order parameter by means of a Hubbard- 
Stratonovich transformation. Contrary to NJL-type 
models, however, the quark self-interactions are not 
free parameters of QCD but generated dynamically and 
driven to criticality by the gauge degrees of freedom. This 
can be understood in simple terms from a renormaliza- 
tion group (RG) analysis of the fixed-point structure of 
four-fermion interactions in gauge theories |20H27j . Such 
a fixed-point analysis also allows for a computation of 
the chiral phase transition temperature as a function of 



the flavor number N{, see Refs. (221 [23]. The results are 
indeed in very good agreement with those from lattice 
simulations [8HTU1 [T3] . Also in agreement with lattice 
studies, the deconfinement and chiral phase transition 
temperature have been found to almost coincide within 
such a first-principles RG setup [H] [TS] . In Ref . [25] it 
has then been shown that the (almost) coincidence of the 
chiral and deconfinement phase transition temperature 
can also be understood by means of an analysis of the 
influence of gluodynamics on the fermionic fixed-point 
structure. 

In particular with respect to the phase boundary of 
QCD at finite temperature and quark chemical poten- 
tial, the interrelation of the chiral and the deconfine- 
ment phase transition is currently under debate, see e. g. 
Ref. [in]- A systematic deformation of QCD represents 
a valuable strategy to gain important insights into this 
question. From the response of the theory to such a de- 
formation, we may then learn something about the un- 
derlying mechanisms at work. For example, varying the 
number of quark flavors, the number of colors, or the cur- 
rent quark masses indicates that the nature of the two 
phase transitions clearly depends on these parameters of 
the theory, see e. g. Refs. [31 l3"Ur[3"5] . More recently, the 
interplay of the chiral and deconfinement phase transition 
has been analyzed by varying the boundary conditions of 
the quark fields in the temporal direction [35]. Such a 
deformation of the theory yields so-called dual observ- 
ables which relate the spectrum of the Dirac operator to 
the order parameter for confinement, namely the dressed 
Polyakov loop QH MM\ ■ 

In this work, we consider a deformation of QCD dif- 
ferent from the ones named above. To be specific, we 
deform QCD by varying the representation of the quark 
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fields. For example, one may consider quark fields in the 
adjoint or in the fundamental representation. From such 
a straightforward deformation it is then possible to gain 
further insights into the mechanisms close to the finite- 
temperature phase boundary of QCD. In fact, basic prop- 
erties of the theory can change when we change the rep- 
resentation of the quark fields: for example, quarks in the 
adjoint representation do not break the underlying center 
symmetry of the gauge sector, whereas the center sym- 
metry is broken explicitly for quarks in the fundamental 
representation. Since center-symmetry breaking is con- 
nected to the question of quark confinement [48] , one may 
expect that a variation of the representation of the quark 
field leaves its imprints in the phase structure of the the- 
ory. In fact, it has been found in lattice simulations of 
SU(iV) gauge theory with adjoint quarks 1!) ">]] that the 
chiral phase transition temperature is significantly larger 
than the deconfinement phase transition temperature, in 
contradistinction to QCD with two (light) quark flavors 
in the fundamental representation. For a study of adjoint 
QCD with Polyakov-loop extended NJL (PNJL) models, 
we refer the reader to Refs. [5!?1 155] , 

Here we pursue the strategy of Ref. [25J and investigate 
the fixed-point structure of quark self-interactions. On 
the one hand, the results of our study provide further 
insights into the interrelation of quark confinement and 
chiral symmetry breaking. On the other hand, our fixed- 
point analysis may help to develop new effective QCD 
low-energy models or to improve existing models. In this 
spirit, the present study also extends previous works in 
which it has been discussed how the low-energy sector 
of QCD can be systematically connected to the QCD 
Lagrangian at high momentum scales within a continuum 
approach [U H3 HJHM \M\ ■ 

This work is organized as follows: In Sect. [IT] we dis- 
cuss general aspects of our field-theoretical setup, with 
an emphasis on the order parameters associated with 
quark confinement and chiral symmetry breaking. The 
fermionic fixed-point structure is then discussed in detail 



in Sect. Ill for general representations of the quark fields. 
In Sect. |IV| we study the partially bosonized version of 
our fermionic ansatz discussed in Sect. |III| The mapping 
between the two formulations is discussed in Sect. IIV Bl 
While the purely fermionic formulation of the matter sec- 
tor already allows us to analyze the mechanisms at work 
at the phase boundary, the partially bosonized version 
allows us to gain access to the hadronic spectrum of the 
theory in a simple manner. In Sect. |IV C[ we then take 
into account l/d(R)-corrections and discuss their effect 
on the finite-temperature phase boundary. Our conclud- 
ing remarks are given in Sect. [V] 



II. GENERAL ASPECTS OF QCD AT LOW 
ENERGIES 



In QCD phenomenology, the quarks are usually as- 
sumed to live in the fundamental representation of the 
underlying SU(A^) gauge group. In the construction of 
general gauge theories, however, the quarks are by no 
means bound to live in the fundamental representation. 
In principle, they may transform according to any irre- 
ducible representation of the SU(^V) gauge group, e. g. 
the fundamental representation (./V-climcnsional) or the 
adjoint representation (TV 2 — 1-dimensional) . On the other 
hand, the gauge degrees of freedom always transform ac- 
cording to the adjoint representation of the gauge group. 

In order to analyze the interplay of the chiral and the 
deconfinement phase transition in QCD, we deform QCD 
by varying the representation of the quark fields and by 
adding a relevant coupling to the theory which can be 
considered as an external deformation parameter. For 
the latter, we choose a four-fermion coupling A^: 



S, 



QCD 



S, 



QCD 



d 4 x X^p (ipCip)' 



where .Sqcd denotes the (classical) action of QCD, see 
e. g. Refs. [55] [5J5] for lattice studies of this class of 
theories. The operator C will be determined below. 
Loosely speaking, such a deformation allows us to "de- 
tune" the chiral and the confining dynamics of the the- 
ory. Due to the additional coupling, the theory ("A^,- 
deformed QCD") now effectively depends on two param- 
eters, namely A^, and AqcdQ In particular, the values of 
chiral low-energy observables, depend on these two pa- 
rameters. Different values of A^, can then be related to 
different values of a given low-energy observable, such as 
the pion decay constant f v . For Aw, = 0, we are left with 
real QCD and the only input parameter is given by Aqcd 
or, equivalently, by the value of the strong coupling a s at 
some (high) momentum scale. In this case, Aqcd solel 
sets the scale for all physical observables O: O ~ Aqcd 

In the following we would like to exploit the depen- 
dence on the two parameters Aqcd and A^, to gain in- 
sights into the relation of the chiral and the deconfine- 
ment phase transition. To this end, we set up a model 
which shares many aspects with the full theory but can 
also be analyzed analytically to a large extent. Following 
Ref. [28] , this allows us to come up with a prediction for a 
phase diagram in the plane spanned by the temperature 
and the pion decay constant f w for different representa- 
tions of the quark fields. This prediction can then be 
tested with the aid of other approaches, such as lattice 
simulations. 



1 Here, we assume that the current quark masses are set to zero. 

2 One could also turn the argument around and, for example, fix 
the scale by choosing a certain value T<j for the deconfinement 
phase transition temperature. The latter then determines the 
scale Aq Cd . 
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We would like to mention that gauge theories with 
an additional relevant parameter, such as a four-fermion 
coupling, have also attracted a lot of attention in re- 
cent years in beyond standard model applications, see 
e. g. Refs. [SU [57]. In particular, the case of SU(2) 
gauge theory with two adjoint quarks is of interest, see 
e. g. Refs. [55H5T] . With regard to our present study, 
a word of caution needs to be added at this point. In 
our numerical analysis in Sects. |Hl| and |IV| we mostly re- 
strict ourselves to the case of SU(2) gauge theory with 
two massless adjoint quarks. We are aware of the fact 
that this theory could already lie in the conformal win- 
dow. In this work, however, we assume that the zero- 
temperature ground state of SU(2) gauge theory with 
two adjoint quarks is governed by dynamical chiral sym- 
metry breaking which would indeed appear to be the case 
within our present approximations. This is in accordance 
with Refs. [57j [62] . As a first step, it is therefore natural 
for us to consider this theory in our numerical studies. 
Even if our present approximations should turn out to 
be insufficient to describe correctly the chiral ground- 
state properties of SU(2) gauge theory with two massless 
adjoint quarks, we still expect that our results will be 
similar for gauge groups of higher rank and broken chi- 
ral symmetry in the zero-temperature limit, see also our 
discussion in Sects. Hill and IIVI 

Before we now study the interplay of the chiral and the 
deconfmement phase transition in detail, we summarize 
a few field-theoretical aspects and explain the general 
setup which underlies our study. 



A. Gauge Sector 

Since we are interested in an analysis of the relation of 
the deconfmement and chiral phase transition in general 
gauge theories (with non-trivial center), such as SU(A) 
and Sp(7V) gauge theories, we need to discuss at least 
some properties of the order parameters for confinement 
and chiral symmetry breaking. In the following we first 
present a few important analytic relations for the con- 
finement order parameter. These relations will play an 
important role in our classification of gauge theories in 
the remainder of this work. 

The deconfmement phase transition in pure SU(iV) 
gauge theories has been studied in great detail. For ex- 
ample, results are available from lattice simulations, see 
e. g. Refs. HHini H3 [32 ESHS1], as well as from func- 
tional continuum methods [33] [70], [71]. A well-known 
order parameter for the deconfmement phase transition 
is the Polyakov loop. The associated Polyakov-loop vari- 
able reads 



1 



L ^ = W) Ve 



j9 So dx o A (x ,x) 



(1) 



notes the bare gauge coupling and V stands for path 
ordering. The Polyakov loop is then given by (trpip). 

Strictly speaking, the Polyakov loop (trpip) is an or- 
der parameter for center symmetry breaking, see e. g. 
Ref. [H]. However, its logarithm can also be viewed as 
half of the free energy F qq of a quark-antiquark pair at 
infinite distance. A center-symmetric confining phase is 
signaled by a vanishing Polyakov loop and implies that 
the free energy of a static fundamental quark (fundamen- 
tal color source) F q ~ (l/2)F qq is infinite. The associated 
quark-antiquark potential is linearly rising for large dis- 
tances in this phase and no string breaking occurs. On 
the other hand, the deconfined phase is associated with 
a finite free energy F q and a finite Polyakov loop, i. e. 
(spontaneously) broken center symmetry, see below. 

Apart from the Polyakov loop, other order parame- 
ters for quark confinement have been introduced, such 
as dual observables [36 . In this work, however, we shall 
mainly consider an order parameter which is closely re- 
lated to the standard Polyakov loop, namely trpI/FKA))]- 
In Polyakov-Landau-DeWitt gauge it has indeed been 
shown that the quantity trFLp[(^o)] serves as an or- 
der parameter for confinement [70] [71]. Here, (Ao) is a 
constant element of the Cartan subalgebra of the gauge 
group and denotes the ground state of the order pa- 
rameter potential in the adjoint algebra, namely the so- 
called Polyakov-loop potential]^] In a one-loop approx- 
imation this order parameter potential has first been 
computed in Refs. (72] [73]. Based on a functional RG 
approach, a non-perturbative study of this potential, in- 
cluding a computation of the phase transition temper- 
atures for several gauge groups, has first been carried 
out in Refs. (33] [7J2 [71]. Since the phase transition 
temperature of a given gauge theory represents a phys- 
ical quantity, it can be easily compared to results from 
other approaches, such as lattice gauge theory, see e. g. 
Refs. [32] [64] [66] [67] . For example, such a comparison 
shows that the RG result for the deconfmement phase 
transition temperature Td for SU(3) is in very good agree- 
ment with results from lattice simulations o 

In the following we shall refer to the temperature as 
the deconfmement phase transition temperature, even if 



where /3 = 1/T is the inverse temperature, e?(F) is the di- 
mension of the fundamental representation of the gauge 
group (e. g. d(F) — N for SU(7V) gauge theories), g de- 



3 Strictly speaking, we have to distinguish between the background 
temporal gauge field and its expectation value (Aq) associated 
with the order parameter for confinement, trpLp[(Ao}]. We skip 
this subtlety here and refer to {Ao) as the position of the ground- 
state of the order-parameter potential if not indicated otherwise. 

4 The deconfmement phase transition has also been studied with 
matrix models, see e. g. Refs. [35 74- 76 . Based on input from 
lattice simulations, on the other hand, ways to improve Polyakov- 
loop potentials widely used in PNJL/PQM-type model studies 
have been discussed in a recent review, see Ref. |77j . For recent 
progress in this direction, we refer to Refs. 1781 1791 . We empha- 
size that our present study relies directly on the order-parameter 
potential spanned by the background temporal gauge field. The 
position (Ao) of the ground state of this potential appears in the 
Feynman diagrams associated with the dynamics in the matter 
sector, see discussion in Sects. Ill B| IIIII and IIVI 
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we study theories with quarks in a representation other 
than the fundamental representation. As we shall briefly 
discuss below, there exist representations for which the 
associated quark-antiquark potential is not linearly rising 
at large distances, independent of the temperature T. 
Below Td, however, center symmetry is restored also in 
these cases. In any case, free (static) color charges are 
always screened to form color-neutral states below Td for 
any representation that is considered in this work. 

Let us now discuss some important properties of the 
quantity 



A 9 So dx Aq(x q ,3) 



(2) 



Here, R denotes the representation of the matter fields, 
e. g. fundamental (F) or adjoint (A), and d(R) is the 
dimension of the representation R. 

Under an arbitrary center transformation of the 
ground state (Ao), 

(A )^(A ) Z , (3) 

the quantity tr^ Lr[(Aq)] transforms as 

tT K L R [(A )] -> z^tr R L R [(A Q )} , (4) 

where A/r denotes the so-called A-ality of the representa- 
tion R. For SU(AT) gauge theories, for example, we have 



z e 



fe 2 ""} 

I J n=0. 



,N-1 



(5) 



Moreover, we have A/r = 1 for R=F. We then obtain 

tr F L F [(A )] -> ztr F L F [(A )] (6) 

and conclude that tr F L F [(Ao)] is an order parame- 
ter for center symmetry breaking and, loosely speak- 
ing, signals confinement of (static) quarks in the fun- 
damental representation, see our discussion above. Note 
that (tr F L F [Ao]) transforms accordingly under arbitrary 
center transformations of the gauge field Aq. Hence, both 
tr F L F [(Ao)] and (tr F L F [Ao]) represent order parameters 
for center symmetry breaking. 

By construction, (Aq) is an element of the Cartan sub- 
algebra. In the following we may therefore parameterize 
it in terms of the generators of this subalgebra: 



Pg(A ) = 2ir 



d(C) 

E 

a=l 



d(C) 

T {a)^{a) =27T^T( a V a )|^|, (7) 
a=l 



where v 2 = 1, the T^'s are the generators of the under- 
lying gauge group in a given representation R. We shall 
refer to the set {4>^} as the coordinates of (A ). 
The center-symmetric phase is signaled by [7T] 

tr F L F [(A )] = (tr F £ P [A ]) =0. 

In the class of Polyakov-DeWitt gauges, it then follows 
that the position (Ao) of the center symmetric ground 
state is uniquely determined by (27J EU [70] 

tr P (i P [(A )n=0, (8) 



where n = 1, . . . , d(C) is the dimension of the associated 
Cartan subalgebra. For SU(A^), we have d(C) = N - 1. 
At (asymptotically) high temperatures, on the other 
hand, we are in the perturbative regime where (Ao) — > 
and tr F L F [(A )] -» 1, see e. g. Refs. [331 [M73]. 
Since tr F _L F [-] is a monotonic function in the domain 
defined by the trajectory of (A ) as a function of the 
temperature T, the quantity tr F L F [(Ao)] is monotonic 
and tr F L F [(Ao)] > in the phase with broken center 
symmetry, see also our discussion below. 

The coordinates {4>^ a '} of the center-symmetric ground 
state, that are determined by Eq. pi), are given by {1/2} 
for SU(2) Yang-Mills theory and {2/3,0} for SU(3), re- 
spectively. From our discussion it is also clear that 
the order-parameter potential in the adjoint algebra is 
periodic^] The lengths of the periods in the various 
directions depend on the eigenvalues of the associated 
generators T^- a K Center transformations of the ground 
state (A ) can now be viewed as discrete rotations of the 
coordinates {4>^} around the center symmetric point. 
For example, we have a reflection symmetry with respect 
to (f> = 1/2 for SU(2). The associated center transforma- 
tion can then be written as follows: 



= 1 



(9) 



with (f) £ [0,1/2]. Under such a center transformation, 
the order parameter transforms according to 

trpiF[(Ao)] = trpLF[0] = cos (7r</>) 

-> tr F L F [(A ) z ] = -tr F L F [(A )] , (10) 

as expected from Eq. Q. Note that we have used Eq. ^ 
to express (A ) in terms of <p. 

For SU(3), center transformations of the ground- 
state (Ao) can be written as rotations by angles of 
27m/3 around the center-symmetric point {2/3, 0}, where 
n = 0, 1,2. Under such transformations, one then finds 
that the order parameter transforms as given in Eq. 

We would like to add that Eq. ^ holds only for the 
center symmetric ground state (A ) for n mod 
For n mod N = and N even, we have 



tr F (L F [(A )] n ) = (-l)^ 



N 1 



(11) 



For odd N and n mod N = 0, on the other hand, we 
have 



tr F (L F [(A )] n ) 



(12) 



Let us now turn to representations other than the 
fundamental one. From Eq. Q, we also observe that 
there may exist representations R of the gauge group for 



5 The order-parameter potential is also invariant under discrete 
rotations about the origin. The corresponding rotation angles 
are determined by the gauge group under consideration. 
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which tr R Z/ R [(^4 )] does not represent an order param- 
eter for center symmetry breaking. To be specific, we 
consider tr R L R [(^4 )] f° r SU(iV) and R=A (adjoint rep- 
resentation). Since we have A/r = (zero iV-ality) in this 
case, we find that tT A L A [(Ao}] transforms as 



tr A L A [(A )} ^ tr A L A [(A ) 



(13) 



under a center transformation. To be more specific, for 
SU(2), we have 

ti A L A [(A )]=tr A L A [cf>} = - [1 + 2 cos (2^0)], (14) 

which is insensitive to arbitrary (center) transformations 
of (A ), see Eq. Thus, tr A L A [(A a )} is not an or- 

der parameter for center symmetry breaking, see also 
Refs. [801 HQ- 

The insensitivity of ti A L A [(A )] and (tr A L A [A }) with 
respect to center transformations is related to the fact 
that quarks in the adjoint representation do not break the 
underlying center symmetry of the gauge group, in con- 
trast to quarks in the fundamental representation. From 
a phenomcnological point of view, there is indeed no 
strict notion of confinement of quarks in the adjoint rep- 
resentation, even in the static limit. In this case, (static) 
quarks can be screened by the gluonic degrees of free- 
dom and form a color-singlet state, as can be seen from 
the decomposition of the tensor product of two adjoint 
multiplets. Therefore a quark-antiquark pair at large dis- 
tances can split up into two singlet states. The associ- 
ated quark-antiquark potential thus flattens at large dis- 
tances and does not rise linearly, as it is the case for 
(static) quarks in the fundamental representation, see 
e. g. Refs. [4^H5TllE0ll8T] for lattice studies. In particular, 
the Polyakov-loop (tr A L A [Ao\) is finite for all tempera- 
tures [80j [81]. Since (tr A L A [A ]) is related to the free 
energy of a static (adjoint) quark, it follows that the free 
energy is finite, even in the center symmetric phase at 
low temperatures. Note that the behavior of the quanti- 
ties (ti A L A [A ]) an d ti A L A [(A }] changes qualitatively 
at T = Td, even though they do not represent order 
parameters for center symmetry breaking. As we shall 
discuss below, this is due to the fact that (ti A L A [A ]) 
and tr A L A [(Ao)] can be related to the order parame- 
ters (tr F .L F L4 ]) an d ti F L F [(A }], respectively. 

Let us close this subsection by summarizing a few use- 
ful relations for the quantity tr R L R [(Ao)]. First, we note 
that the order parameter tr F L F [(Ao)] can be related to 
the standard Polyakov-loop via the Jensen inequality. 
For a given concave function /(•), we have /((■)) > (/('))• 
For example, this yields 



tv F L F [(A a )} > (tr F T F L4 ] 



(15) 



for SU(2) and SU(3) gauge theory in the deconfincd 
phase [331 [7D1 EH- We emphasize that this inequality 
does not hold for general gauge groups and representa- 
tions R, since it requires that tr R T R [-] is a concave func- 
tion in the relevant domain. Provided that (A ) (T) lies 
sufficiently close to the origin (e. g. for sufficiently large 



temperatures T), however, the inequality may hold for 
any gauge group and representationp] 



In addition to Eq. ( 15 ), we have the following two sim- 
ple but useful inequalities: 

\tr R (L R l(A )}^\ < jj— , (16) 



< 



1 



<^e[tr R (T R [<A ) 



l )l < 1 (17) 

n ~ (d(R))™ 1 ' 



d(R) 1 nA VLN "~ (d{R)) n 

which follows from the generalized triangle inequality, 
and 

1 

for neN. 

Finally, we evaluate the quantity tr R L R [(Ao)] for 
specific configurations of the ground-state (Ao). At 
very high temperatures T ^> Td, we have (Aq) — > 
and tr R L R [(^4 )] L independent of the gauge group 
and the representation R. In the low-temperature phase 
(T < Td), however, the value of tr R L R [(Ao)] depends 
on the gauge group and the representation R. For exam- 
ple, we have ti F L F [(A }] = for SU(iV) gauge theories 
for T < Td. For the adjoint representation, on the other 
hand, we find 

tv A (L A [(A )] n ) ' ' 



(d(A)) n (N 2 - 1)« 
with n mod N ^ 0. For n mod N = 0, we have 



tv A (L A [(A )r) = 



(d(A)y 



(18) 



(19) 



In SU(iV) gauge theories, the relation (18) follows 



straightforwardly from the fact that the tensor product 
of the triplet and the anti-triplet can be decomposed into 
the adjoint multiplet and a singlet, TVigiA = (jV 2 -l)® 1. 
Using that the character of the product representation is 
given by the product of the characters of the representa- 
tions, we find 

d(A) tr A (L A [<i4o)]) = \d(F) tr F (L F [(A )])| 2 -l , (20) 

and similar relations for ti A (L A [(Ao)] n ) with n > 1 (n E 
N). 

Up to this point, we have discussed that below Td 
the quantity tr R L R [(A )] is zero for R=F but negative 
for R=A. Depending on the representation R, however, 
it can also assume positive values in the center symmetric 
phase. For example, let us consider the ten-dimensional 
representation (R = 10) of SU(3). We then fincQ 

d(10)tri (£io[<4))]) 

= [d(F)tr F (T F [(A )])] 3 

-2|d(F)tr F (T F [(A )])| 2 + l = l (21) 



6 In general, a subdomain around the origin can be found such that 
trnLn[-] is a concave function. For a more detailed discussion 
of the relation | |15| between the order parameters trp£F[(^lo}] 
and (trpLp[ylo]}, we refer the reader to Ref. 128) . 

7 Here, we have used that the tensor product of the triplet and the 
sextet can be decomposed into a singlet and decuplet. Moreover, 
the tensor product of two triplets can be decomposed into a 
sextet and an anti-triplet. 
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for T < Td, and similar relations for tri (.i/io[(A))] n ) 
with n > 1 (ra G N). Note that the 7V-ality of 
the ten-dimensional representation of SU(3) is zero as 
well, A/io = 0. We add that also higher-dimensional rep- 
resentations exists for which we have tr R LR ,[(Ao}] = 
for T < Td, even if the iV-ality of the representation is 
zero. 

Depending on the representation R, we have seen that 
the quantity tr R L R [(^4 )] can be zero, positive or nega- 
tive in the center symmetric phase, i. e. for T < T^. This 
allows us to classify gauge theories. In Sect. we shall 
see that this classification is to some extent related to 
the question whether the chiral phase transition temper- 
ature is larger or smaller than the deconfinement phase 
transition in a given gauge theory 



B. Matter Sector 

Let us now discuss our field-theoretical setup in the 
matter sector. Up to this point, our statements concern- 
ing the pure gauge sector are exact and can be obtained 
analytically. They only rely on the basic assumption that 
we work in the class of Polyakov-DeWitt gauges. For our 
analysis of the matter sector, we employ the following 
ansatz for the quantum effective action T: 

r[^,(A )]= f d 4 x{z^j(i?) + g l0 (A })i; 



[(#) 2 - (^r 75 V>) 2 ] } , (22) 



where is the wave-function renormalization of the 
quark fields. In the present work, we restrict ourselves 
to N{ = 2 massless quark flavors with d(R) colors. The 
t,-'s represent the Pauli matrices and couple the spinors 
in flavor space. 



Our ansatz ( 22 ) for the matter sector is perturbatively 



non-renormalizable, as it is the case for the NJL model. 
Therefore we define it with an UV cutoff A which then 
represents an additional parameter of the model. This 
setup also implies that the regularization scheme belongs 
to the definition of the model. The role of A for the fixed- 
point structure will be discussed in detail in Sect. |IV| 

For our study of the RG flow of the four-fermion 
coupling we employ the so-called Wetterich equa- 
tion |82j . The latter is an RG equation for the quan- 
tum effective action. In this approach, the effective ac- 
tion r depends on the RG scale k (infrared cutoff scale) 
which determines the RG 'time' t = ln(fc/A) with A be- 
ing a UV cutoff scale. For reviews on and introductions 
to this functional RG approach, we refer the reader to 
Refs. HZHHH1]. 

Concerning the background field (A)), we will 
not make use of the approximation tr R L R [(/lo)] = 
(tr R L R L4 ]) which underlies most PNJL/PQM model 



studiesQseee. g. R e fs. [3 H [5J [7J HQ [HQ [S3 El [MM] . 
Although this assumption is convenient and opens up the 
possibility to incorporate lattice results for (tr R L R L4 ]), 
it may be problematic in quark representations other 
than the fundamental one and away from the limit of 
infinitely many colors, d(R) — > oo, see our discussion 
below. For our analytic studies, we shall rather make 
use of the exact relations given in Sect. |II A| For our 
numerical evaluation of the quantum effective action, 
we then use the numerical results for (Aq) from a non- 
perturbative first-principles RG study of the associated 
order parameter potential in Polyakov-Landau-DeWitt 
gauge, see Refs. [33 [TD]. Since it has been found for fun- 
damental matter that PNJL/PQM- type model studies 
are sensitive to different parameterizations of the poten- 
tial for (trpiF^o]): see e. g. Ref. [98;, it is in fact impor- 
tant to analyze at least some of the consequences arising 
from the approximation tr R L R [(^4o)] = (tr R L R [Ao]) un- 
derlying these model studies. 

In general, our ansatz (22) in the matter sector can be 



considered as the leading order in a systematic deriva- 
tive expansion. The associated expansion parameter is 
the anomalous dimension r)^ — —df In Z^p of the quark 
fields. This "parameter" is small as has been found in 
various previous studies [201 EH (SSI HDD]. In fact, it is 
identical to zero when we consider the four-fermion cou- 
pling in the so-called point- like limit, A^,(|p| <C fc), see 
e. g. Ref. [37]. This is true even if we had allowed for dy- 
namical gauge degrees, provided that one considers the 
class of Landau gauges |101) , such as Polyakov-Landau- 
DeWitt gauge. As we have discussed above, the latter 



gauge is implicitly assumed in our work, see Sect. II A 



Apart from an expansion in derivatives, the effective 
action can be expanded in operators, such as n-fermion 
operators. Regarding four-fermion operators, we note 
that our ansatz (22) for the effective action is not com- 



plete with respect to Fierz transformations even in the 
limits (A ) ->■ and T -> 0, see e. g. Refs. [SI]- 
1351 I3T1 UOlj . For example, we have dropped a so-called 
axial-vector channel interaction which would also con- 
tribute to the RG flow of our coupling Aj, associated with 
a scalar-pseudoscalar channel. From a consideration of 
a Fierz-complete basis, however, we only expect quan- 
titative corrections to our results presented here. The 
main qualitative aspects are expected to persist since 
the general structure of the loop integrals remains un- 
changed [35]. For a Fierz-complete study of RG flow of 
four-fermion couplings in QCD, we refer the reader to 
Refs. [3TH3J [SSI H3 EH]- Regarding the role of higher 
fermion operators, e. g. 8-fcrmion operators, we note that 
it can be shown that these operators do not contribute to 
the RG flow of the four-fermion couplings in the point- 
like limit [ST]. Beyond the point-like limit, however, these 



8 Note that there are also PNJL/PQM- type model studies which 
do not use this approximation but consider an integration over 



the group SXJ(N), see e. g. Refs. |61|93|. 
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higher-order operators may very well contribute to the 
flow of the four-fermion interactions, see, e. g., our dis- 
cussion in Sect. IV and Ref. [2"7j . 

In the subsequent section we will show that the purely 
fermionic formulation of our ansatz ( 22 ) for the matter 



sector is convenient for a general discussion of the inter- 
play of the chiral and the deconfinement phase transi- 
tion, independent of the fermion representation. In or- 
der to compute low-energy observables, however, a purely 
fermionic formulation may not be the first choice since 
this requires to resolve the momentum-dependence of the 
fermionic vertices. In this partially bosonized 

formulation of our ansatz ( 22 ) might be better suited. 



Such a formulation of the effective action can be obtained 
straightforwardly from a Hubbard-Stratonovich transfor- 
mation of the underlying path integral and yields 

r fc 0,V,*,<4)>]= / d 4 x{z^(i$ + g l0 (A ))i> (23) 

}■ 



1 



m 2 $ 2 



where the auxiliary bosonic fields <f> T — (a, ir) mediate 
the interaction between the fermions. Here, we consider 
these bosons to be composites of fermions which do not 
carry any internal color or flavor charges: a ~ (ipip) an d 
7? ~ {jpf^^ip). Chiral symmetry breaking is now signaled 
by a non- vanishing expectation value of the a held. Since 
the mass parameter fa 2 describes the curvature of the 
chiral order parameter potential at the origin, the sign 
of fa 2 is related to the question whether chiral symmetry 
is broken in the ground state or not. For fa 2 < 0, we 
necessarily have (<j) ^ 0. 

In the following we choose the initial conditions for the 
various couplings in Eq. ( 23 ) such that 



lim m 2 > . 



lim Zq, = , 

fc->A 



lim Zjh = 1 . 



Together with the identity 



A - 12 
va 



(24) 



the ansatz ( 23 1 can then be mapped onto the ansatz ( 22 ) 



at the initial UV scale A. Thus, only the ratio of the 
Yukawa coupling h and the mass parameter fa acquires a 
physical meaning. In particular, we observe that a large 
(i. e. diverging) four-fermion coupling signals the onset 
of chiral symmetry breaking, since it can be related to 
a change in the sign of the parameter fa 2 . In fact, the 
two criteria are equivalent in the large-d(R) limit due 
to the absence of fluctuation effects of the Goldstone 
modes [37]. In any case, we conclude that the fixed- 
point structure of the coupling A,/, (or, equivalently, of 
the couplings fa and h) is directly linked to the ques- 
tion of chiral symmetry breaking in the IR limit. In the 
following we analyze how this fixed-point structure is re- 
lated to the order parameter for center-symmetry break- 
ing, namely trpLp[(ylo)]. This will eventually allow us to 
gain insights into the relation of chiral symmetry break- 
ing and center-symmetry breaking at finite temperature. 



T > k (A ) 




<A ) > 



Figure 1. Sketch of the f}\^ -function of the four-fermion 
coupling A^, for vanishing temperature (black/solid line), fi- 
nite temperature and {Ao) = (red/dashed line), and fi- 
nite tem per ature and (Ao) > (blue/dashed-dotted line), 
see Eq. ( |25[ ). The arrows indicate the direction of the RG 
flow towards the infrared. The figure has been taken from 
Ref. [27]. 



III. DYNAMICAL LOCKING MECHANISM 
AND THE FERMIONIC FIXED-POINT 
STRUCTURE 

Let us begin with an analysis of the interplay of center 
symmetry breaking and chiral symmetry breaking using 
the purely fermionic formulation of the effective action, 
see Eq. (22). Our discussion follows closely the analysis 
in Ref. 28]7 where the deformation of the fermionic fixed- 
point structure due to the presence of confining dynamics 
has been analyzed for fundamental quarks in detail. 

In our study we consider the value of the background 
field (Ao) as an external input which is given by the 
ground state of the corresponding order parameter poten- 
tial. As discussed above, the position (Ao) of the ground 
state is then directly related to our order parameter for 
center-symmetry breaking, namely trp£F[(^o)]- Along 
the lines of Ref. [55] , the RG flow equation of the four- 
fermion coupling in the point-like approximation can be 
computed for quarks living in any representation R. We 
find 



Px =d t X^ = (2 + 27fy)A,/, 



(25) 



d(R) 



-^+^)£<i r, M,,rww, 

where the dimensionless renormalized coupling A^ is de- 
fined as 



A^, — Z^ 2 k 2 \ijj . 



(26) 



For example, we have d(R = F) = N for quarks in the 
fundamental representation and d(K = A) = N 2 — I for 
quarks in the adjoint representation. For convenience, 
we have introduced the eigenvalues vi of the hermitian 
matrix given in Eq. ([7]): 



vf 1 =spec{(T<V% | v 2 = 1} 



(27) 



The coupling depends on the background field (Aq) 
and the dimensionless temperature r = T/k. The thresh- 

(F) 

old function l\ describes a regularized one-particle 
irreducible (1PI) Feynman diagram with two internal 



fermion lines. The definition of this function can be found 
in, e.g., Ref. El 



The RG flow equation ( [25] ) has two fixed-points]^] a 
GauBian fixed-point (A^ = 0) and a non-trivial fixed- 
point X^(t, (Aq)), see Fig. [TJ The non-Gaufiian fixed- 
point can also be computed analytically]^] 



a;(t, (a )) 



/ 1 

A;(0,0) 1 + — ^(-d(R))" [tr R (L R [(^o>] n )+tr R (4[(A )]")] (l + ") 

\ ^ ' n—1 



(28) 



where 



a;^a;(o,o) = 



6n 2 
(2d(R) + 1) 



(29) 



the scale for all chiral low-energy observables O: 

A! 



f 



0(\7 - A 



■0/ ' 



(30) 



where do is the canonical mass dimension of the observ- 
able O and the critical exponent O is defined as 



e 



d\,i 



2. 



(31) 



Thus, we have (i(R)A^ const, for d(R) ->• oo. Note 
that we have dropped terms depending on 77^ on the 
right-hand side of Eq. ( [28] ) . As discussed above, this is 
not an approximation in the point-like limit. 

Before we turn to the case of finite temperature, we 
briefly discuss a few basic aspects of the zero-temperature 

limit. In order to solve the RG flow equation §25§, we For details, we refer the reader to Ref. [27| . In the fol- 
have to choose an initial value A^ v at the scale k = A for lowm S we fix A ^ V > A; at T = 0. The value of A^v 
the coupling A,/,. For A^ v < A^, we find that the four- 
fermion coupling approaches the Gaufiian fixed-point in 
the IR limit, i. e. the theory becomes non-interacting and 
chiral symmetry remains intact. For Aj^ v > A^, on the 
other hand, we observe that the four-fermion coupling 
grows rapidly and diverges at a finite scale fcsB- This 
scale signals the onset of chiral symmetry breaking. Be- 
low this scale, the point-like approximation is no longer 
justified: the formation of a quark condensate and the 
appearance of Nambu-Goldstone modes requires that we 
resolve the momentum dependence of the four-fermion 
coupling in this regime. For example, this can be done 
by means of partial bosonization techniques, see Sect. |IV| 
In any case, the chiral symmetry breaking scale /csb sets 



then determines the symmetry breaking scale fcsB and, in 
turn, the values of the chiral low-energy observables. For 
our study of finite-temperature effects and effects from 
the confining dynamics parameterized by the background 
field (j4o), w e leave our choice for A^J V unchanged. This 
ensures comparability of our results at zero and finite 
temperature for a theory defined by a given value of A^J V . 

At finite temperature and finite (Aq), the fixed-point 
structure of the theory is deformed compared to the zero- 
temperature limit. For illustration purposes, we begin 
with a brief discussion of the case with vanishing gluonic 
background field. In this case, the pseudo fixed-point is 
shifted to larger values at finite temperature, A^,(r, 0) > 



XI, see Eq. j28|. For a given initial value A^ v > A^, , this 



9 At finite temperature, the non-Gaufiian fixed-point is rather a 
pseudo fixed-point, i. e. the fixed-point inherits an implicit scale 
dependence from the dimensionless temperature r = T/k as well 
as a dependence on {Ao). Moreover, the line of pseudo fixed- 
points A^(t, {Aq)) does not represent a separatrix in the (A^,,r)- 
plane. However, it represents a strict upper bound for the sepa- 
ratrix in this plane I27| . 
10 The fixed-point value is not a universal quantity as can be seen 
from its dependence on the regularization scheme. However, the 
statement about the existence of the fixed-point and its qualita- 
tive dependence on the temperature and (Aq) is universal. 



implies that a critical temperature T x exists above which 
chiral symmetry is restored, see Ref. [27] for a detailed 
discussion. Strictly speaking, the critical temperature T x 
is defined to be the temperature for which 1/A^, — > 
for k — > 0. From the RG flow equation ( [25] ) , one then 
obtains a simple analytic expression for T x : 



T x = 



A 



1 - 



K) , (32) 



which is accordance with our general statement in 
Eq. ( 30 1 . To derive this expression, we have assumed 
that T/A< 1. 
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Let us now turn to the case of finite (Aq). For fermions 
in the fundamental representation, for example, we have 
tr F (L F [(A )] n ) —> in the center symmetric phase 
for n e N and d(F) = N -> oo, see Eq. p]). Thus, the 
temperature-dependent corrections to A^,(r, (Ao)) van- 
ish identically and we have A^(r, (A )) = A^(0, 0) for 
T < Tj. We shall refer to this as a locking mecha- 
nism for the chiral phase transition [55] ■ Loosely speak- 
ing, this means that the chiral phase transition is locked 
in due to the confining dynamics. For T > Td, we 
have trFLp[(A )) > and the fixed-point is again shifted 
to larger values. As pointed out in Ref. [55J, this implies 
that T x > T d in the limit N -s- oo, see also Refs. [T1UD2]. 
In the case of adjoint fermions and c?(R) 3> 1, the 
temperature-dependent corrections in Eq. ( |28[ ) do not 
vanish identically on all RG scales k for T < Td since 
we now have trA(TA[(Ao)]") < for these tempera- 
tures, see Eq. (18). Therefore the (global) sign of the 



temperature-dependent corrections changes compared to 
the case with (Ao) — 0. This implies that the pseudo 
fixed-point is shifted to smaller values at finite tempera- 
ture rather than to larger valuesj^j] For fixed T and k — > 
(i. e. r — > oo), the pseudo fixed-point approaches 



X*Jt^oc,(A )) 



A* 



d(A) 



(33) 



A^ from be- 



For d(A) ->• oo, we have X*Jt -> oo, (A )) 
low. Since trA-kA[(Ao}] -> 1 for T > T d , the fixed-point 
is "released" and shifted to larger values. For a given ini- 
tial value A^ v > A^,, it then follows again that T x > Td 
for d(A) — > oo. 

This analysis can in principle be repeated for any 
fermion representation, including fermion representa- 
tions for which tr R L R [(A )] > for T < T d , such as the 
ten-dimensional representation, see Eq. (21). In the lat- 



ter case, the temperature-dependent corrections do not 
vanish for large values of d(R). However, the finite- 
temperature shift of the fixed-point is still suppressed by 
a factor of l/d(R) compared to the case with (Ao) = 0. 
As a consequence, the chiral phase transition tempera- 
ture is increased, but it is not necessarily pushed above 
the deconfinement phase transition temperature. There- 
fore a strict statement about the relation of the chiral 
and the deconfinement phase transition cannot be made 
for this class of theories, not even in the large-d(R) limit. 

From Eq. ( 25 ) , we can derive an implicit equation 
for the chiral phase transition temperature T x pro- 
vided that we neglect a possible RG scale dependence 
of tr R i R [(A )]: 



rp2 _ 1 

x ~~ r R (T x ) 



A 



\UV 



0(a7-a;),(34) 



11 Note that an external magnetic field deforms the fixed-point 
structure in a similar way 103, 104 . 



where 



d(R)7r 2 ^ n 

x ' 71 — 1 



oo 1 

£-(-d(R)r[tr R (T R [(A )n 



tr R (L R [(A )n] . (35) 



In order to derive this equation, we have again assumed 
that T/A < 1. For (A ) -> 0, we have V R {T) = 1, 
as expected from Eq. (32). Since tr R L R [(A )] depends 
on Td, Eq. (34) relates the chiral phase transition tem- 



perature T x to the deconfinement phase transition tem- 
perature Ti- 

For fermions in the fundamental representation 
and N -> oo, we have V F (T) = for T < T d 
and V F (T) > for T > T d . Thus, no finite solution 
of Eq. ( 34 1 exists for T x < T d . In accordance with our 
fixed-point analysis, we conclude that T x > Td. 

For fermions in the adjoint representation and d(A) 3> 
1, we have T A (T) « -l/d(A) for T < T d . As in the case 
of fundamental matter fields, this implies again that T x > 
Td. For fermion representations with < 'Pr(T) < 1 
for T < T d , we observe that the chiral phase transition 
temperature is still shifted to larger values compared to 
the case with (Ao) = 0. However, a definite statement 
about the temperature-order of the two phase transitions 
cannot be made in this case. 

Let us now turn to the case of finite values of d(R). 
For fermions in the fundamental representation, we then 
find -Pf(T) = 1/N 2 > for T < T d . For T > T d , 
"Pf(T) increases monotonically from V-piT^) = 1/N 2 to 
Vf(T — > oo) — > 1. In fact, right above the deconfinement 
phase transition, the quantity Vp(T) increases rapidly 
since trpTp[(Ao)] increases rapidly. Since 'Pf(T) is finite 
for all temperatures, we find a regime where T x < Ti 
for Aj^ v /A; > 1. In the left panel of Fig. [2J we show 
our numerical results for T x /T d as a function of A^| v /A^. 
For larger values of A^ v /A^,, a window in parameter 
space opens up in which the chiral and the deconfinement 
phase transition (almost) coincide. In the limit N — > 00, 
this locking window extends down to A^ v /A,0 = 1, as 
illustrated by a comparison of our results for N = 2 
and TV = 3 in Fig. [2] Note that the locking window 
for A^ v /A.0 can be related to a locking window for low- 
energy observables, such as the pion decay constant. Wc 
shall come back to this in Sect. IIVI 

In case of adjoint fermions and finite d(A), we have 
Va(T) < for T < T d . For example, we have V A (T) = 
-1 for N = 2. For T > T d , V A (T) increases rapidly 
and changes its sign. For T >• Td, it then ap- 
proaches V A (T) = 1. Since we have V A (T) < even for 
finite N, we find that the chiral phase transition temper- 
ature is larger than the deconfinement phase transition 
temperature, independent of our choice for A^ V /AJ, > 1 
and N > 2, see right panel of Fig. [2] for our numerical 
results for TV = 2. Note that this observation is com- 
patible with lattice results of SU(2) gauge theory with 
two adjoint quarks, see Refs. [4T)rl51| . In our analysis, it 
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N = 2 (fund.) 
N = 3 (fund.) 



1.4 



1.6 1.8 



2.2 2.4 



N = 2 (adj.) 
N = 2 (fund.; 



1.4 



1.6 1. 

k UV, 



2.2 2.4 



Figure 2. Left panel: Phase diagram in the plane spanned by the temperature and the rescaled coupling A^ v /AiJ, for Nf = 2 
massless quark flavors in the fundamental representation and N — 2 colors (red/solid line) as well as for N = 3 colors 
(blue/dashed line), see also Ref. [2S]. Note that there is no splitting of the phase boundary (i. e. T x ~ Td) for small A^ v in the 
large- iV limit, see Eq. (34 1 and discussion thereof. Right panel: T x /Td as a function of A^ v /Ay, for Nf = 2 massless quarks in 
the fundamental representation (N = 2) (red/solid line) as well as for quarks in the adjoint representation (blue/dashed line). 



can be traced back to the deformation of the fermionic 
fixed-point structure in the presence of gauge dynamics. 

To obtain the numerical results in Fig. [2j we have em- 
ployed data for (Aq)(T) as obtained from an RG study 
of the associated order parameter potential for SU(2) 
and SU(3) Yang-Mills theory [33l [70]. However, we 
did not take into account the back-reaction of the mat- 
ter fields on the order parameter potential associated 
with (Aq). In the case of fundamental matter, we expect 
that this back-reaction will shrink the size of the locking 
window since it further increases the quantity Vf(T) at 
low temperatures. For adjoint quarks, the back-reaction 
will also increase Va(T). Nevertheless, it may remain 
negative over a wide range of temperatures. Therefore 
we may still have T x > T d for all values of A^ v /A^ > 1, 
at least for N = 2. 

Let us add a word of caution on the treatment of the 
quantity tr R L R [(A )) in standard PNJL/PQM model ap- 
proaches. In these studies, one relies on the assumption 
that tr R L R [(A )} = (tr R L R [A }}. For (tr R L R [A ]), one 
then uses lattice data as input. Whereas such an ap- 
proach would lead to similar conclusions for fundamental 
quarks ((trpLF[A ]) > and trF-kp^o)] — 0)> tne 
uation is different for adjoint quarks. In the latter case, 
we have (trA-^AlA)]) > but trA-kA[(^4o)] can assume 
both positive and negative values as discussed above. 

Before we now enter the discussion of the RG flows of 
the partially bosonized formulation of the matter sector, 
we would like to comment on the number of parameters 
in our study. Up to this point, our discussion suggests 
that our study only relies on a single parameter in the 
matter sector apart from the UV cutoff A, namely on the 
initial value A^ v . Strictly speaking, however, the non- 
trivial fixed-point of the four-fermion interaction is an 
artifact of our point-like approximation. With the aid 



of the partially bosonized formulation, we will resolve 
part of the momentum dependence of the four-fermion 
interaction. We will then find that the matter sector 
depends on three parameters: the Yukawa coupling h, 
the bosonic mass parameter m and the UV cutoff A, 
see Eq. (23 1. This is a substantial difference to, e. g., 



fermion models in d < 4 space-time dimensions, where 
we only have a single parameter in both formulations, see 
e. g. Ref. [105 . There, the non-trivial fixed-point of the 
four-fermion coupling can be mapped onto a correspond- 
ing non-trivial fixed-point in the plane spanned by the 
renormalized Yukawa coupling h and the dimensionless 
renormalized bosonic mass parameter m. In our case, the 
role of the non-trivial fixed-point in the purely fermionic 
formulation is taken over by a separatrix in the (h 2 , m 2 )- 
plane in the partially bosonized formulation. The shift 
of the non-trivial fixed-point of the four-fermion coupling 
due to the gauge dynamics then turns into a correspond- 
ing shift of this separatrix. The mapping between the two 
formulations is discussed in detail in the subsequent sec- 
tion. Being aware of this subtlety, the discussion of the 
fermionic fixed-point structure is still useful and nicely 
illustrates the mechanism underlying the interplay of the 
chiral and the deconfinement phase transition. 



IV. PARTIAL BOSONIZATION AND THE 
LARGE-d(R) EXPANSION 

A. Gap Equation 

In this subsection, we briefly discuss how our study 
of fermionic RG flows is related to the gap equation for 
the fermion mass in the large-d(R) limit. For related 
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QCD reviews on Dyson-Schwinger equations, we refer the 
reader to Refs. [17)BriT09] . 

Starting from the partially bosonized action given in 
Eq. ( 23 ) , we can derive the gap equation for the vac- 
uum expectation value $o = (( cr )i0) of the scalar fields 
and the fermion mass, respectively. To this end, we first 
consider the so-called classical action S ~ which 
appears in the functional integral. Since the fermions 
appear only as bilinears in the action S, these fields can 
be integrated out straightforwardly. From the resulting 
expression, we obtain the (fully) bosonized effective ac- 
tion Tefcr, 7r], which is a highly non-local object. From 
the stationary condition, 



ST 



B 0-, 7T j 



8a 



= 0, 



(36) 



we then find the gap equation for (a) P] 

A) J (2-) 3 

-G|; i) (A 2 ,.9(Ao),(a))l0(A 2 -p 2 ), (37) 



G^(p 2 ,g{A ),{a}) 



where 



G^>(p',g(A ),(a)) = 



(38) 



(u n +g(A Q )) 2 +p 2 + h 2 (a) 2 
and v n = (2n + \)nT. The trace Tr is defined as follows 



Tr-- - =tr R T ^ 



(39) 



n— — oo 



In Eq. (37) we have dropped the trivial solution (a) = 0. 



The integral on the right-hand side of the gap equation 
represents a Feynman integral with two internal fermion 
lines and two (Yukawa) vertices. In order to regularize 
this integral, we have employed a regularization schemj^] 
which corresponds to the one used to derive the RG flow 
equation for the four-fcrmion coupling A^, in the previous 
section, see Eq. (25). The structure of the loop integral 



in the gap equation and on the right-hand side of the 
flow equation (25) is indeed identical^] However, the 
prefactor on the right-hand side of the gap equation is 
only correct in leading order in the large-d(R) expan- 
sion, in contrast to the associated prefactor in the flow 
equation (25) of the four-fermion coupling Our gen- 
eral arguments concerning the relation of the chiral and 
the dcconfincmcnt phase transition in the previous sec- 
tion are not affected by this prefactor. The latter plays 



12 Here and in the following we assume that the ground state is 
homogeneous. 

13 Often, a sharp cutoff is used to regularize the gap equation . In 
this case, the A-dependent term in the square brackets in Eq. | |37| l 
is absent, see Ref. |27| . 

14 Recall also that hj^/fh^ can be identified with Aj^ v = A^, a- 

15 Loosely speaking, the trace tr^ yields a factor of d(R). 



only a qualitative, but no quantitative role. Therefore 
our findings concerning the interplay of the chiral and 
the deconfinemcnt phase transition can be obtained from 
the gap equation (37) as well, as it should be. In fact, 



the fixed-points of the (dimensionless) four-fermion cou- 
pling can be viewed as critical values for the dimension- 
less quantity A 2 fo 2 /m 2 . This follows also from our dis- 
cussion below Eq. ( 23 1 . We refrain here from discussing 
this further. For a detailed discussion of fermionic RG 
flows and their relation to the gap equation, we refer the 
reader to, e. g., Ref. [37]. 



B. RG Flow at Large d(R) 

Let us now discuss the fixed-point structure and the 
locking mechanism in the partially bosonized formula- 
tion of the matter sector, see Eq. (23). This formulation 



has the advantage that it allows us to systematically re- 
solve the momentum dependence of the four-fermion in- 
teraction by means of a derivative expansion, see e. g. 
Ref. [27]. Eventually, this allows us to relate the initial 
value A^J V of the four-fermion coupling to physical low- 
energy observables, such as meson masses and the pion 
decay constant f n . The phase diagrams in the (T, A^ v ) 
plane can then be translated into phase diagrams in, e. g., 
the (T, /„.) plane. In other words, the partially bosonized 
formulation gives us access to the phase with broken chi- 
ral symmetry in the ground state. 

From the effective action ( p3| , we obtain the RG flow 
equations for the partially bosonized formulation. In 
leading order of an expansion in powers of d(R), we find 
the following equations for the chirally symmetric regime: 



d(R) 



1=1 

m< = o , 

d t h 2 = (2»ty + m)h 2 , 



d(R) 



(40) 

(41) 
(42) 



d t m 2 = ( m -2)m 2 + ^Y, li(r,0,^\)h 2 , (43) 



<9tA$ = 2?7$A^ 



i=i 



d(R) 

1=1 



(44) 



where 77$ 



2£2 



-d t \nZ^, h 



= k- 



! Z^m 2 



and A* = Z^ 2 A$. The 



threshold functions can be found in App.[A|and Ref. [24j . 
We add that we do not distinguish between wave-function 
renormalizations longitudinal (zl, z\) and transversal 
(Z^, Z$) to the heat-bath. In the following we iden- 
tify the corresponding wave-function renormalizations, 
Z\ = Z^ = Z^ and Z\ = Zi = Z*. In Ref. [ED], 
it has indeed been found for the case (An) = that the 
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difference is small at low temperatures and only yields 
mild corrections to, e. g., the thermal masses close to 
and above the chiral phase transition. In our study of 
the partially bosonized formulation, we also include the 
RG flow of the four-boson coupling A$ associated with an 
additional term ~ (A$/8)$ 4 in our ansatz (23). Since we 
have $ ~ (tjjOtp), this type of interaction parametrizes 
higher-order fermionic self-interaction terms. These in- 
teractions are generated dynamically in the RG flow due 
to Yukawa-type quark-meson interactions. The initial 
value of the associated coupling is set to zero in our stud- 
ies, i. e. A$ = at k — A. This allows us to map the par- 
tially bosonized theory onto our purely fermionic ansatz 
for the matter sector at the initial RG scale A. 

In the large-d(R) limit, the flow of the four-boson cou- 
pling (and also of higher bosonic self-interactions ~ $ 2ra ) 
does not contribute to the RG flows of Z$, Z^, h and m, 
at least in the chirally symmetric regime. This corre- 
sponds to the fact that the RG flow of the four-fermion 
coupling is decoupled from the RG flow of fermionic n- 
point functions of higher order, such as 8-fermion inter- 
actions. 

Since we consider the large-(i(R) limit in this subsec- 
tion, we only have purely fermionic loops appearing on 
the right-hand side of the flow equations. We would like 
to stress that the large-<i(R) expansion should not be con- 
fused with the widely used local potential approximation 
(LPA) where the running of the wave-function renormal- 
izations is not taken into account. 

Let us now relate the RG flow of the partially 
bosonized formulation to the RG flow of the purely 
fermionic formulation. Using Eq. ( 24 ) together with the 



flow equations (42) and (43), we recover the flow equa- 
tion ( p5| ) of the four-fermion coupling X^> in the large- 
d(R) limit, i. e. 



= d f X 



d(R)->oo 



(45) 



d(R)->oo 



see also Eq. (63). Thus, only our choice for the ra- 
tio h 2 /m 2 at the initial RG scale determines whether 
chiral symmetry is broken in the IR limit (fc — > 0). Since 
the flow equation for the ratio h 2 /m 2 is identical to the 
one for A^, we already anticipate that our statements 
concerning the temperature-order of the chiral and the 
deconfinement phase transition still hold in the partially 
bosonized formulation, see also Ref. [28 . 

After having shown the equivalence of the RG flow 
of h 2 /m 2 and A^, in the chirally symmetric regime, we 
now discuss the number of parameters in the matter sec- 
tor of our model. Relation (45) seems to suggest that 
we only have one parameter, namely the ratio h 2 /m 2 
at the initial RG scale. Indeed, the value of the sym- 
metry breaking scale /csb depends only on our choice 
for h 2 /m 2 at the initial RG scale. This suggests that 
a non-trivial IR repulsive fixed-point also exists in the 
plane spanned by the couplings h 2 and m 2 . From the 
above set of flow equations, however, we read off that 
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Figure 3. RG flow for fundamental fermions and N — 3 in 



leading order in the l/d(_R)-expansion in the (ft 



-piano 



(at zero temperature). The red (solid) line represents the sep- 
aratrix (critical manifold). The arrows indicate the direction 
of the RG flow towards the infrared, see text for an interpre- 
tation. 



the system has only a Gaufiian (non-interacting) fixed- 
point (h 2 ,m 2 )Gauii = (0,0), but no non-Gaufiian fixed- 
point. This seems to be in contradiction to Eq. (45) 



and to our results from the purely fermionic formulation. 
Apart from the Gaufiian fixed-point, we also observe that 
a separatrix exists in the plane spanned by h 2 and m 2 . 
The latter separates the (h 2 , ro 2 )-plane into two disjunct 
regimes, see Fig. [3] In the large-d(R) limit, the func- 
tional form of the separatrix can be computed analyti- 
cally. At T = 0, we find 



scp 



(m 2 ) 



3vr 2 m 2 
d(R) 



= A 



ip,oo ' 



(46) 



where A,0 ;OO is the value of the fixed point A5, in the 
large-d(R) limit. Choosing initial conditions (h\,m\) 
in the domain to the left of the separatrix, we find that 
the system flows into the regime with m 2 < 0, in which 
chiral symmetry is broken in the ground state. On the 
other hand, the system remains in the chirally symmetric 
regime, if we initialize the flow in the domain to the right 
of the separatrix, see Fig. [3] Loosely speaking, we have 
found that the separatrix takes over the role of the non- 
Gaufiian fixed-point A^ which is present in the point-like 
approximation of the purely fermionic formulation. 

To further clarify the fate of the seemingly missing non- 
trivial fixed-point in the (/i 2 , m 2 )-plane, we briefly con- 
sider the case 2 < d < 4. In this case, a non-trivial fixed- 
point indeed exists in NJL-type and Gross-Neveu-type 
models, see e. g. Ref. 105 . This follows immediately 
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k a ~*Zl ZI h 



2^2. 



from a consideration of the RG flow of the dimensionless 
renormalized Yukawa coupling h 2 

(d- 



h- 



d t h 2 



77$ + 2r]^)h 



(47) 



This differential equation has a Gaufiian fixed-point and 
a non-Gaufiian fixed-point h 2 for 2 < d < 4 since 77$ ~ 
h 2 and 77$ > 0, see Refs. |105j . In the (ft 2 , m 2 )-plane, 
we therefore have a non-trivial fixed-point with an IR 
attractive and IR repulsive direction for 2 < d < 4. This 
non-trivial fixed-point represents the intersection point of 
two separatrices in the (h 2 , m 2 )-plane and corresponds to 
the non-trivial fixed-point of the associated four-fermion 
coupling. For d — > 4, this fixed-point then merges with 
the Gaufiian fixed-point. 

The non-existence of the non-Gaufiian fixed-point 
in d = 4 implies that the Yukawa coupling h and the UV 
cutoff A should be considered as parameters of the the- 
ory, in addition to the ratio h\/m\. In fact, for any finite 
UV cutoff, we can still define a critical value for the ratio 
h 2 /m 2 ~ A.0 above which spontaneous chiral symmetry 
breaking occurs in the long-range limit (fc — > 0). How- 
ever, since no non-trivial fixed-point with an IR attrac- 
tive direction exists in the (h 2 , m 2 )-plane, the value of the 
Yukawa coupling at the symmetry breaking scale &sb de- 
pends (strongly) on the initial conditions for the bosonic 
mass parameter and the Yukawa coupling itself. Once the 
system enters the regime with m 2 < 0, the RG flow of the 
Yukawa coupling effectively "freezes" . In this low-energy 
regime, the fermions acquire a mass and fermionic loops 
are therefore generically suppressed, see also discussion 
below. Hence we have three parameters in our simplified 
model ansatz for the matter sector, namely h\/m\, h\ 
and A. We would like to emphasize that there is indeed 
only a single parameter in 2 < d < 4, as discussed in 
detail for the Gross-Neveu model in Ref. |105j . 

Let us now analyze the dynamics at finite T and (Aq). 
Our discussion of Eq. (45) and of the RG flow in the 
(/1 2 , m 2 )-plane at zero temperature already suggests that 
our general arguments concerning the relation of the de- 
confinement and the chiral phase transition in Sect. Ill 
are still valid. This is not too surprising: the point- 
like approximation in the purely fermionic formulation 
is a reasonable approximation in the chirally symmetric 
regime where the bosonic mass parameter m 2 is large 
over a wide range of scales and therefore suppresses the 
non-trivial momentum-dependence of the vertices^] In 
any case, we now have to study the behavior of the sepa- 
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Figure 4. Sketch of the RG flow in leading order in the 



l/d(R)-expansion in the (h 



-plane. The black dot de- 



notes the Gaufiian fixed-point. The separatrices are sketched 
for three different cases: vanishing temperature (black line), 
finite temperature and (.Ao) = (red/dashed line), and finite 
temperature and (Ao) > (blue/dashed-dotted line). The 
dependence of the separatrices on the temperature and (Ao) 
reflects the behavior of the non-Gaufiian fixed-point of the 
four-fermion coupling, see Fig. [I] The arrows to the left and 
to the right of the separatrices indicate the direction of the 
RG flow towards the infrared, respectively. 



one-dimensional manifold as it is the case at zero temper- 
ature. It rather represents a two-dimensional manifold. 
Again, the functional form of this critical manifold can 
be computed analytically. For r = T/A <C 1, we find 



cp.( m V) 



1 - 7T 2 V R {T)T 2 ' 



(48) 



where T*r,(T) is defined in Eq. (35). We observe that 



the shape of the critical manifold depends on the tem- 
perature and the order parameter for center symmetry 
breaking, see Fig. [4] 

The critical manifold allows us to define a necessary 
condition for chiral symmetry breaking at finite temper- 
ature. Solving Eq. (48) for r, we obtain r scp .(m 2 , h 2 ). 

scp. for a given set of initial val- 
the theory necessarily approaches the 



Choosing now r < 7, 
ues (h 2 



A ' 



nn. 



regime with broken chiral symmetry in the IR limit. 



For r > T m 



on the other hand, the theory remains 



ratrix in the [h , to )-plane for finite temperature T and in the chirally symmetric regime. For a given value of 



finite (Ao). To this end, we may even consider the di- 
mensionless temperature t = T/k as an additional cou- 
pling of the theory. Thus, the separatrix is no longer a 



the UV cutoff A and (hj^, mjA, the quantity r scp . is there- 
fore nothing but the dimensionless chiral phase transition 
temperature, r scp . = T x /A. In fact, T x = At scp . agrees 
with the result from Eq. (34). Thus, our general state- 



16 Recall that the bosons mediate the interaction between the 
fermions in the partially bosonized formulation. In this spirit, 
the boson propagators parametrize the momentum dependence 
of the four-fermion coupling, see e. g. Ref. [27| for a detailed 
discussion. 



ments in Sect. |III| concerning the interplay of the chiral 
and the deconfinement phase transition still hold. 

Let us now discuss how our phase diagrams in 
the (T, A^ v )-plane can be translated into phase diagrams 
in, e. g., the (T, / 7r )-plane. To this end, we need to follow 
the RG flow down to the long-range limit. As discussed 
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A [MeV] U [MeV] 

Figure 5. In the left panel, we show the phase diagram for two massless fundamental quarks and N = 3 in the plane spanned 
by the rescaled temperature T x /Td and the value of the pion decay constant j n at T — 0. In the right panel, the corresponding 
phase diagram for two massless quark flavors in the adjoint representation and N = 2 is shown. In both panels, the results 
from the large-d(R) approximation are given by the red (solid) line, whereas the blue (dashed) line depicts the results from our 
study including corrections beyond the large-d(R) limit. 



in Sect. II B the mass parameter m 2 assumes negative 
values in the regime with broken chiral symmetry in the 
ground state and the vacuum expectation value (<&) = $o 
becomes finite. It is therefore convenient to study the RG 
flow of $o an d A$ rather than that of m 2 and A$. The 
flow equation of 3>o can be obtained from the stationary 
condition: 



d 
dt 



d 



1 



2^2 



-m $ 



1 



= 0. 



(49) 



To be specific, we find the following RG flow equations 
for the regime with broken chiral symmetry in the ground 
state: 

d(R) 



d t h 2 



3tt 2 ^ 
i=i 



0, 



m)h 2 , 

-2)4? 

s d(R) 



(r, m 2 , 



(51) 
(52) 

(53) 



d t U = 27?*A* - -^Yl 'a (r.m^vj^Dfc 4 , (54) 



i=i 



where $q 



k 2 Z$$ 2 and the (dimensionless) renormal- 
ized constituent quark mass reads 



m\ = hH 2 



In the following we will identify the pion decay con- 

1 /2 — 

stant /,r with $o- The (dimensionless) renormalized 
meson masses are given by 



and 



Since we are working in the large-d(R) limit in this sec- 
tion, the latter do not appear explicitly on the right side 
of the flow equations. 

Recall that the scale for m q and m a is set by the 
symmetry breaking scale fcgB which is set by our choice 
for h\/m 2 A . The role of the Yukawa coupling (as an ad- 
ditional parameter) becomes now apparent from the re- 
lation 



= A*$ 



which follows from the flow equations of the couplings. 
Since the flow of the Yukawa coupling is not governed 
by the presence of a non-trivial IR attractive fixed-point, 
its value depends on A;sb and the initial value h\, as 



(50) discussed above. Therefore the ratio m 2 /m 2 depends on 



our choice for h\. On the other hand, the initial value 
of the coupling A$ does not represent a free parameter 
of the theory. It is set to zero at k = A and therefore 



generated dynamically in the RG flow, see also Eq. (23 1 



Using the flow equations (40 1- (44) and (50 1- (54 1, we 



can now proceed and compute the phase diagram in the 
plane spanned by the temperature and the value of the 
pion decay constant at T = 0. In Fig. [5] (left panel) we 
show our results for quarks in the fundamental represen- 
tation and N = 3. For adjoint matter and N = 2, our 
results can be found in the right panel of Fig. [5j To ob- 
tain these results, we have used A = 1 GeV. Moreover, 
we have again employed the data for the ground-state 
values of (^4o) a s obtained from a RG study of SU(A r ) 
Yang-Mills theories [531170] , 

In the case of fundamental matter and N = 3, we ob- 
serve that the upper end of the locking window (Td ps T x ) 
roughly coincides with the physical value of the pion 
decay constant, provided that we fix the initial condi- 
tion of the Yukawa coupling such that m q sa 300 MeV 
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for fa « 90 MeV, see left panel of Fig. [5| This obser- 
vation is in accordance with results from lattice simu- 
lations and general expectations. For / w < 30 MeV, we 
find T x < T(j. More precisely, we observe that T x ~ / w for 
small values of f„. For / w > 100 MeV (m q > 350 MeV), 
we then have T x > In this regime, the quarks are 
very heavy and the two phase transitions are disentan- 
gled. Concerning the role of the Yukawa coupling, we 
find that the lower end of the locking window is shifted 
to smaller values of / w when we increase the initial value 
of the Yukawa coupling. Moreover, we find that the size 
of the window does not strongly depend on our choice 
for h\. This is not unexpected since we have found in 
our analysis of the fermionic fixed-point structure that 
the size of the locking window is solely related to the 



d t h 2 = {2 n4 , + m )h 2 



2 1 
7? d(R) 



d(R) 



J2lff\r,0,^\<l>\,m 2 )h\ (56) 



i=i 



d t m 2 = { m -2)m 2 - —h{T,m 2 )\ 



2tt 2 



d(R) 



1=1 
3 

<9 t A$ = 2r?$A$ + -^1 2 (t, m 2 )\% 
d(R) 



^4 F) (r,O,^|0l)^ 4 , 



(57) 



(58) 



i=i 



\UV 



However, the transla- 



as 



value of the ratio h\/m 2 

tion of the upper and lower end of the locking window in 
A^ v -space into values of physical observables does indeed 
depend on our choice for both h\ as well as h\/m 
discussed above. 

For adjoint matter and N 
find that T x > Td, even for very small values of /„■. We 
refer to Fig. [5] for our results for N = 2. To obtain 
these results, we have used h\ = 3. However, T x > Td 
holds for arbitrary values of in the large d(R) limit, as 
suggested by our fermionic fixed-point analysis. In fact, 
our results in the large-d(R) limit are in accordance with 
our results in Fig. [2j as it should be. For increasing f n , 
we observe that the chiral phase transition temperature 
increases further. Thus, we have T x > Td for all values 
of U. 

Finally we would like to add that it is also possible 
to tune the parameters h\/m\ and 



In regime with broken chiral symmetry ($ 7^ 0), the 
flow of the couplings is determined by the following equa- 
tions EH 



d(K) 



2 as well as N = 3, we V<s> = tt^ E •^S.C 7 "' m v ^H)^ 2 > 



(59) 



i=i 

,2 _ /o„ i „ \u2 



d t h z = (2^+r? 4 )/i 2 



1 1 

^ d(R) 



d(R) 



E [ 3 4 F l B) ( T ' m q'^l^l' m ^) 
1=1 



(60) 



such that we 



obtain T x /Ta ss 7.8 for N — 3, as found in lattice simu- 
lations [50] of adjoint QCD without A^-deformation. Of 
course, this requires that the UV cutoff A is adjusted to 
larger values in order to ensure that T/A is sufficiently 
small for the temperature range under consideration, see 
Ref. [52] for a PNJL model study in a mean-field approx- 
imation. 



C. RG Flow Beyond the Large-d(R) Approximation 

In the following we study the robustness of our re- 
sults of the previous sections with respect to l/d(R)- 
corrections. This includes an analysis of the role of 
Goldstone-mode fluctuations which are absent in the 
large-G?(R) limit. 

Our RG approach allows us to systematically in- 
clude l/c?(R)-corrections. Due to the one-loop structure 
of the Wetterich equation, these corrections correspond 
to 1PI diagrams with at least one internal boson line. In 
the chirally symmetric regime ($ = 0), we then find the 
following set of equations: 



d t <P 2 = -( m + 2)#§ + _i_ Zl ( r>m 2) + _^_ Zl ( T)TO 2) 

8 d(R) h 2 

1=1 A * 

9 3 
<9 t A$ = 2t7$A$ + ^Z 2 (t,to^)A| + — 1 2 (t, ml)\% 



famish*. (62) 



i=i 

The threshold functions^] can be found in App. [X] and 
Ref. [23]. For simplicity, we do not include the run- 
ning of the fermionic wave-function renormalization in 
the present study, although it can be taken into account 
straightforwardly, as illustrated in, e. g., Refs. [20 ] [24 ] [99 ] 
1100] for the case (Aq) = 0. As discussed in the previous 
subsection, this is not an approximation in the large-d(R) 



2 
3^ 



d(R) 

i=i 



(55) 



17 In the flow equations for the Yukawa coupling and the bosonic 
wave-function renormalization, we have dropped terms propor- 
tional to $o- Concerning the Yukawa coupling, it has been 
found that these terms only yield mild (quantitative) correc- 
tions 20, 110, 111 . With regard to the bosonic wave-function 
renormalization, these terms are of crucial importance for an ac- 
curate computation of the critical exponents 1991 11121 which is 
beyond the scope of the present work. 

18 Note that the functions and depend implicitly 
on 
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limit. Beyond the large-d(R) limit, it has been found in 
Refs. [20j [24j |99l 1100) that the anomalous dimension 77^, 
is still small. This can be traced back to the fact that the 
running of Z<$, is solely governed by 1PI diagrams with at 
least one internal boson and fermion line. Such diagrams 
are parametrically suppressed in the regime with broken 
chiral symmetry due to the large mass of the fermions, 
but they are also suppressed in the chirally symmetric 
regime due to the large mass of the bosons. As a con- 
sequence, the running of Z^ only yields mild corrections 
to the symmetry breaking scale &sb- ^ n the following we 
will only take into account l/c?(R)-corrections in those 
RG equations which are also non-zero in the large-d(R) 
limit. The inclusion of the running of is left to future 
work. 

Using the flow equations of the Yukawa coupling and 
the bosonic mass parameter in the chirally symmetric 
regime, we can study again the RG flow of the ra- 
tio h 2 /m 2 . We now find 



<2(R) 



2?r 



3 ( h 2 



1=1 



, d W 

— — Vz (FB) f7 



,m 2 ) 



{=1 



(63) 



Using = (h 2 /to 2 ), the first term on the right-hand side 
as well as the terms in the second line can be straightfor- 
wardly identified with terms appearing in the RG equa- 
tion of the four-fermion coupling A,/,, see Eq. (25 1. These 



are the leading order terms of the large- g?(R) expan- 
sion. The second term on the right-hand side corresponds 
to a l/c?(R)-correction and effectively couples the flow 
of h r/m (~ four-fermion coupling) to the flow of the 
four-boson coupling (~ 8-fcrmion coupling) . Since it can 
be shown that the RG flow of fermionic self-interactions 
is fully decoupled in the point-like limit [27], this term 
resolves (part of) the momentum dependence of the four- 
fermion interaction. The expression in the third line also 
represents a 1 /<i(R)-correction and can be traced back to 
the running of the Yukawa coupling. Without the terms 
in the third line, it is not possible to reproduce the prefac- 



tor of the term ~ in the flow equation ( 25 ) in the limit 



771 ^> 1 (point-like limit). As pointed out in Ref. [28) . 
this can be seen immediately from the following relation 

1 

— 2 "i v , ..- q) 



2 (m»l) 



ID- 



4 P) (r,m 2 vm 



For finite m 2 , the expression in the third line on the right- 
hand side of Eq. ( 63 1 also resolves part of the momentum 



structure of the four-fermion vertex beyond the point-like 
limit. 

Let us now discuss our results for the phase diagrams in 
the (T, / w )-plane beyond the large-d(R) limit. In Fig. [5] 
we show our results for quarks in the fundamental repre- 
sentation and N = 3 as well as for quarks in the adjoint 



representation and N = 2. We have chosen these repre- 
sentations and values for N since they play a prominent 
role from a phenomenological point of view. For quarks 
in the fundamental representation, we observe that our 
results agree with those from our large- iV study, at least 
on a qualitative levelj^] This means we still have three 
distinct regimes: one regime with T x < Td for small val- 
ues of /„■, one regime with T x m Td (locking window), 
and a regime with T x > Ti for large values of f„. Also, 
the size of the locking window is roughly the same as in 
the large- N approximation. However, the lower and the 
upper end of the window have been shifted to larger val- 
ues of f„. The locking window begins at / T ps 100 Me V 
and ends at f v ~ 150 MeV. Thus, the physical value of 
the pion decay constant is slightly below the lower end 
of the locking window. 

For quarks in the adjoint representation and N = 2, 
we find that our results are less strongly affected by cor- 
rections arising beyond the large-d(A) approximation, 
see Fig. [5] (right panel). To be specific, we observe 
that T x > T d for f„ > 0, even if we take l/d(R)- 
corrections into account. The results only differ with 
respect to the slope of the chiral phase transition tem- 
perature as a function of f v . As in the case of fermions 
in the fundamental representation, the slope is steeper 
in the large-<i(A) limit. We conclude that fluctuations of 
the Nambu-Goldstone modes tend to lower the sensitivity 
of T x on /„.. 

The results for adjoint quarks in Fig. [5] have been ob- 
tained by choosing h\ — 3 for the initial value of the 
Yukawa coupling. The value of the pion decay constant 
can then be varied by varying only the initial value of 
the bosonic mass parameter m\. As in the case of funda- 
mental quarks, it is in principle possible to fix the initial 
condition for the Yukawa coupling by requiring that the 
constituent quark mass assumes a given value for a given 
value of the pion decay constant. For adjoint quarks, 
we refrain from fixing the initial condition h\ in this way 
but rather illustrate how our results depend on the choice 
for h\, see Fig. [6j We observe that the dependence of T x 
on becomes stronger for larger values of h\. Most 
importantly, however, we find that T x > Ti for N = 2, 
independent of our choice for Ha > 0. We stress that the 
mechanism underlying this observation is the deforma- 
tion of the (fermionic) fixed-point structure due to the 
presence of the confining gauge dynamics. 

Let us finally comment on the order of the chiral phase 
transition in the (T, f w ) phase diagram. In Ref. [28], it 
was found for fundamental fermions and N = 3 that the 



19 Note that we have fixed the initial value of the Yukawa cou- 
pling by requiring that m q ft! 300 MeV for f w ft 90 MeV. The 
same initial value for the Yukawa coupling has then been used 
to compute the phase transition temperature for all other values 
of f„ as well. Thus, we have only varied the initial value of the 
bosonic mass parameter m 2 to change the value of f^. Recall 
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Figure 6. Ratio T x /Td of the chiral and the deconfine- 
ment phase transition temperature as a function of the zero- 
temperature value of the pion decay constant f n for two mass- 
less adjoint quarks and N — 2. The various lines illustrate the 
dependence of our results on the initial condition (UV value) 
for the Yukawa coupling. The results have been obtained 
for ft a = 2, 3, 4 (from bottom to top). 



chiral phase transition is of first order within the lock- 
ing window. To be more precise, we observe that the 
chiral phase transition is of first order for 100 MeV < 
fa < 150 MeV for N = 3. Above and below the lock- 
ing window, the chiral phase transition is of second or- 
der. In particular, the observation of a first-order region 
might be a shortcoming of our approximations: we have 
simply used the data for (Aq)(T) from a study of pure 
SU(3) Yang-Mills theory, but neglected the back-reaction 
of the matter sector on the confinement order parameter. 
Within the locking window, the first-order phase transi- 
tion in the gauge sector induces a first-order chiral phase 
transition. As argued in Ref. [28], a first-order chiral 
transition may still occur in the (T, / 7r )-plane, even if 
we go beyond the present approximation. However, this 
would then require that the confinement order parame- 
ter rises rapidly for T > T^. A test of this conjecture is 
beyond the scope of the present work and left to future 
studies. 

For adjoint matter and N — 2, we observe that T x > 
Id for all values of > 0. Therefore the dynamics at 
the chiral phase transition is less affected by the confining 
dynamics. Loosely speaking, the latter only pushes T x 
above Within the present approximation, we there- 
fore find that the chiral phase transition is of second or- 
der for all values of > 0. This result is consistent with 
lattice simulations for TV = 2, see Refs. fl9l [ST)] . 



V. CONCLUSIONS AND OUTLOOK 

In the present paper we have analyzed the interplay 
of the chiral and the deconfinement phase transition in 



gauge theories with matter fields in different representa- 
tions, with an emphasis on quarks in the fundamental 
and the adjoint representation. To this end, we have 
computed phase diagrams in the plane spanned by the 
temperature and the pion decay constant using a simple 
ansatz for the quantum effective action. This ansatz al- 
lowed us to study the fixed-point structure in the matter 
sector analytically. In particular, it opened up the pos- 
sibility to analyze the impact of the confinement order 
parameter on the chiral fixed-point structure. The lat- 
ter is directly related to the order parameter for chiral 
symmetry breaking. 

For theories with quark fields living in a given rep- 
resentation R, we have found that the interplay of the 
chiral and the deconfinement phase transition clearly de- 
pends on the sign of the quantity trRiR^Ao)] in the cen- 
ter symmetric phase. The relation of trRLpi[(Ao)] to the 
standard Polyakov loop (for a given representation R) 
has been discussed in detail in Sect. |II A] To be specific, 
our fixed-point analysis suggests that T x > for adjoint 
quarks, at least in the large-d(R) limit. This observation 
is in accordance with results from lattice simulations [4"9T - 
I51j . For quarks in the fundamental representation, our 
findings are also compatible with lattice QCD studies [5J- 
fTO] H3] , first-principles continuum studies [12l [15] 1113] , 
and earlier analytic (model) studies ]T] [55] 1102) . 

We have also investigated how robust our predic- 
tions for the (T,f % ) phase diagram are, once l/d(R)- 
corrections are taken into account in the matter sector. 
Such corrections are associated with fluctuations of the 
Nambu-Goldstone modes of the theory. For quarks in 
the fundamental representation, we have found that the 
locking window (T x wlj) is shifted to larger values of 
but remains finite. At the physical point (/„• ~ 90 MeV), 
we have T x < T<j. For adjoint quarks and N = 2, we 
have found that T x > for all values of f n > 0, even 
if l/d(R)-corrections are taken into account. In this re- 
spect, the finite-temperature dynamics of gauge theories 
with adjoint matter appear to be distinct from gauge 
theories with fundamental matter, at least for TV = 2. 
Recall that this observation is also consistent with lat- 
tice studies of SU(2) gauge theory with two flavors of 
adjoint quarks [4"9"r[ST] . 

We would like to add that our fixed-point analysis 
can also help to guide the development of QCD low- 
energy models in the future. To be specific, we have 
used the order-parameter potential spanned by the back- 
ground temporal gauge field. Moreover, we have not em- 
ployed the assumption trpLR [(A )] = (trpjLpifAo]) which 
is often used in PNJL/PQM-type model studies. In- 
stead, we have considered the quantity trr<LR[(^4o)] m 
the present study. This corresponds to working in a 
specific class of gauges, namely the class of Polyakov- 
DeWitt gauges. Our study suggests that the assump- 
tion trR,Ljj,[(Ao)] = (trR,Lri[Ao]) is justified in a mean- 
field approximation (N — > oo). For finite V (and, in par- 
ticular, for quarks in representations other than the fun- 
damental one), however, the situation may change. De- 
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pending on the representation R, the sign of trR,LR[(^4o)] 
can be different below and above the phase transition, 
whereas (tr R L R L4 ]) can be defined to be positive for all 
temperatures. 

Of course, the present analysis can be improved in 
many ways. For example, one may consider to take 
into account the back-reaction of the matter fields 
on the quantity tr R Lp. [(Aq)]. Such contributions will 
push tr R L R [(v4o)] to larger values in the low-temperature 
phase. We expect that this will weaken the mechanisms 
governing the dynamics in our present study. For fun- 
damental quarks, for example, this may shrink the lock- 
ing window. For adjoint quarks, on the other hand, the 
quantity tr R L R [(Ao)] may still be negative over a wide 
range of temperatures. Therefore, T x > may persist 
for N — 2, even if we take these back-reactions into ac- 
count. In any case, we have presented first predictions for 
the (T, /tt) phase diagram and our analysis reveals a sim- 
ple mechanism governing the interplay of the chiral and 
the deconfinement phase transition. It would be interest- 
ing to see whether and how this mechanism persists in 
the presence of a finite quark chemical potential and/or 
a finite external magnetic field. The latter deforms the 
fermionic fixed-point structure in a way [1031 1104] which 
is indeed reminiscent of the deformation discussed here 
for adjoint quarks. 
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specific, we have chosen 
'k 2 



-1 9(k z -p z )=p'r B 



IT 

for the bosons, whereas we have chosen 




k 2 



,(A1) 



(A2) 



for the fermionic degrees of freedom. 

Now we define the threshold function Ai\ j_ . This func- 
tion represents a 1PI diagram with two internal fcrmion 
lines and contributes to the RG flow of the bosonic wave- 
function renormalization. Those threshold functions, 
which are not defined in this appendix, can be found 
in Refs. (2U E2 [28] ■ 

(F) 

To define the threshold function A4\ j_, it is convenient 
to introduce a dimcnsionless propagator for the fcrmions: 

1 



G^(x 0l oj) 



x Q + x(l + r^) 2 + uj 



(A3) 



where x — p 2 /k 2 . In terms of this propagator, the 
threshold function entering the anomalous dimension of 
the bosons can be written as follows: 



OO 

= (d-l)r £ 



2x 

d-l \dx 2 
d + 1 / d 



dxx 2 dt I x(l+r^,)G^,(xg , cu) x 
I 

(1 + r^)G^(xQ , lj) 



d — 1 \ dx 
^ 2x ( d 2 



1 V dx- 



G^(x$,oj) 



—G^,{xi,uj) 



(A4) 



Appendix A: Threshold functions 



In the computation of the RG flow equations, a reg- 
ulator function needs to be specified which determines 
the regularization scheme [82] . Here, we have used a lin- 
ear spatial regulator function for the bosonic as well as 
for the fermionic degrees of freedom |114H118] , To be 



where Xq = (y n + 2ttt^l) 2 and v n — (2n + 1)ttt. Here, 
the derivative dt with respect to the regulator function 
is defined as follows: 



(A5) 



We do not display terms cx 770 since we have not taken 
into account these contributions in our numerical analy- 
sis. 
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